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ABSTRACT 

>- ■ We present new 2.2 /xm diffraction-limited images from the W. M. Keck 10 



0^ 



m and Gemini 8 m telescopes of the cool Galactic Center sources, IRS IW, 5, 8, 



' low, and 21 along with new proper motions for IRS IW, lOW and 21. These ob- 

^ servations were carried out to test the bowshock hypothesis presented by Tanner 

^ ■ et al. as an alternative to a very recent (10"^ yr) epoch of star formation within the 

■ tidal stream of gas and dust known as the Northern Arm. Resolved asymmetric 
p i' structure is detected in all the sources, with bowshock morphologies associated 
Q ■ with IRS IW, 5, 8 and lOW. For IRS IW and lOW, there is an agreement be- 
^ ■ tween the position angle of the asymmetry and that of the relative velocity vector 
^ ■ of the near-infrared source with respect to the Northern Arm gas strengthening 

the bowshock hypothesis. We therefore conclude that the observed morphology 
^ . is indeed a bowshock generated by sources plowing through the Northern Arm. 

■ Furthermore, the large extent of the resolved structures (310-1340 AU) along 
with their luminosities (^lO'*"^ Lq) suggests that their central sources are Wolf- 
Rayet stars, comparable to the broad He emission-line stars, which have strong 
winds on the order of 1000 km s~^. The bowshock morphology, along with the 
proper motion measurements, provide three-dimensional orbital solutions for this 
enigmatic class of objects; IRS IW and lOW have orbital planes that are con- 
sistent with that of the putative clockwise plane which has been proposed as a 
solution for the He I emission-line stars. While these observations eliminate the 
need to invoke star formation within the Northern Arm, they increase by 14% 
the total known population of massive, young stars with strong winds, whose 
origin remains unexplained in the context of the nearby supermassive black hole. 
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Subject headings: Galaxy: center — infrared: stars 

1. Introduction 

Within the Galaxy's central cluster there are a number of enigmatic sources (IRS IW, 
2, 5, 8, low and 21) that have eluded classification for almost three decades (Rieke & Low 
1973, Becklin and Neugebauer 1975). They are all spatially coincident with the Northern 
Arm, a tidal stream of dust and gas that has been mapped at both near-infrared and radio 
wavelengths and appears to be infalling and/or orbiting around the Galaxy's central super- 
massive black hole (Lo & Claussen 1980; Scrabyn & Lacy 1985; Lacy et al. 1979; Herbst et 
al. 1993; Paumard, Maillard & Morris 2004). These sources also all have spectral energy dis- 
tributions that peak near 10 fim, high luminosities (L~10^ Lq, Becklin & Neugebauer 1978), 
and, in the cases of IRS IW and 21, nearly featureless near-infrared spectra (Moultaka et al. 
2004; Krabbe et al. 1995) and polarization properties that cannot be accounted for by the 
intervening ISM (Eckart et al. 1995; Ott et al. 1999). Based on these observed properties, 
many have suggested that these Northern Arm sources are very recently formed stars that 
are being transported into the central core along with the Northern Arm inflow (Krabbe et 
al. 1995; Blum et al. 1996; Ott et al. 1999). As pointed out by Tanner et al. (2002), this 
scenario is problematic since the short Northern Arm timescales require that the stars would 
have to have formed, or to have begun forming, prior to the infall of the Northern Arm. In 
this case one would not expect all of the stars to still be embedded in the gas of the Northern 
Arm, given that the newborn stars will immediately follow ballistic orbits, while the gas in 
which they formed will be subject to strong, additional, non-gravitational forces. 

We recently argued, based on a simple radiative transfer model of the spatially resolved 
2-25 nm images and photometry of IRS 21, that all the Northern Arm sources are massive 
stars experiencing rapid mass loss, as either Wolf-Rayet (WR) or asymptotic giant branch 
(AGB) stars (Tanner et al. 2002). As the Northern Arm flows past these stars, the interac- 
tion between the material in the Northern Arm and the intrinsic stellar winds generates a 
bowshock which produces the observed spatially resolved emission (Ott et al. 1999; Tanner 
et al. 2002). While the resolved structure around IRS 21 lacks the characteristic horseshoe 
morphology, this could simply be a projection effect, if the relative velocity between the 
outflow source and the Northern Arm is oriented near the line of sight. One appeal of the 
bowshock model is that it explains the association of these sources and the Northern Arm 
without requiring that the stars form within the Northern Arm before or during the infall 
toward the central black hole. A bowshock model instead allows for the much more plausible 
circumstance that the sources are simply in the path of the infalling Northern Arm gas and 
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dust, allowing them to be members of the observed population of stars seen in and around 
the central parsec. This hypothesis is supported by recent high-angular-resolution images of 
IRS 8, which show a clear bowshock morphology around this star (Geballe et al. 2004). 

In this paper, we present new near-infrared observations for the remaining Northern 
Arm sources, IRS IW, 5, and lOW, as well as a new analysis of IRS 21 and 8 in order to 

test the bowshock hypothesis. Section 2 describes the data sets that were obtained for this 
work, Section 3 reviews the techniques used in the image construction and deconvolution, 
and Section 4 presents the measured morphologies and proper motions for the surveyed 
sources. Section 5 makes the case that these sources arc Wolf-Rayet stars - young, massive, 
post-main sequence stars - generating bowshocks as they plunge through the Northern Arm, 
and it discusses the implications of their source classification. 

2. Observations 

2.1. Speckle Imaging with Keck I/NIRC 

Speckle imaging observations of IRS IW, 5, lOW, and 21 were obtained in the K band- 
pass (Ao=2.2 //m, AA=0.43//m) using the W. M. Keck I 10-meter telescope and the facility 
near-infrared camera (NIRC; Matthews & Soifer 1994; Matthews et al. 1996). The camera 
has a pixel scale of 20.400±0.042 mas/pixcl (Ghcz ct al. 2004) and a correspondingly small 
field of view (FOV) of 5"x5", limiting observations to one target at a time. Each object 
was observed in sets of ~200x0.14 second exposures. The long-exposure seeing averaged 0'.'6 
during these observations. 

All sources except IRS 5 were observed on multiple nights (see Table 1) as part of a 
proper motion campaign (Ghez et al. 1998, 2000, 2004). Measurements of the 5"x5" field 
centered on Sgr A* began in 1995, providing a eight-year time baseline for measurements of 
the proper motion of IRS 21. Beginning in 1998, additional images were obtained to mosaic 
the central stellar cluster frames together with those that contained the SiO masers, IRS 7 
and lOEE, in order to register the radio and infrared reference frames (see Ghez et al. 2004), 
resulting in a four-year time baseline for IRS IW and lOW. 

2.2. AO Imaging with Gemini/Hokupa'a-|-Quirc^ 

We have also utilized public Gemini Adaptive Optics (AO) 2.2 /xm images^ of IRS 
8 collected on 2000 July 02 using the Hokupa'a AO system and the QUIRC near-infrared 
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camera with a FOV of 20" x20" and a plate scale of 19.745±0.090 mas/pixel. The plate scale 
estimate for these images is based on the average and standard deviation of the independently 
determined plate scales estimated using the NIRC plate scale and the five offsets from IRS 
16C to IRS 16NW, 16SW, 16NE, 33N, and IRRl, compared to those predicted by Ghez 
et al. (2004). In these observations, USNO 0600-28577051, which has an R (0.7/im)-band 
magnitude of 13.8 and is 19" away from Sgr A*, served as the natural guide star. 

3. Data Analysis 

3.1. Image Construction 

The final images for all the observations are produced with the same basic two steps: 
1) each frame is sky-subtracted, flat-fielded, bad-pixel-corrected, and resampled by a factor 
of two (except for the Gemini data), and 2) the images at a given telescope pointing are 
combined by shift-and-adding (SAA) on the brightest source located near the center of the 
field of view. For the speckle images, the centroid of the dominant speckle within the speckle 
interferometry pattern is chosen and for the adaptive optics images, the centroid of the 
star over a region corresponding to the diffraction limit was used. These SAA maps have 
Strehl ratios of ~2-5% and 15% for the Keck speckle images and the Gemini AO images, 
respectively, and it is these images that are deconvolved (see Section 3.2) to study the 
morphology of the Northern Arm sources. 

For the proper motion analysis, the Keck speckle SAA maps at different telescope point- 
ings arc assembled using all the bright sources in the overlap regions between pointings 
(typically 20-30 sources that are brighter than K~14 mag). The final mosaics cover roughly 
10" X 10" . Source positions, both final positions and those used to register the different point- 
ings, are based on peaks in the cross-correlation of the SAA maps with the diffraction-limited 
PSF core as described by Ghez et al. (1998); here, a cross-correlation threshhold of 0.7 is 
used. The positions of all the identified sources within the mosaiced images covering eight 
separate epochs are aligned with respect to a reference image (1999 July) by minimizing 
the net offset of all the stars within the mosaic. The NIRC data set provides 21, 10, and 7 
independent astrometric position measurements for IRS 21, IW, and lOW, respectively. 



^ Based on observations obtained at the Gemini Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the 
Gemini partnership: the National Science Foundation (United States), the Particle Physics and Astronomy 

Research Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the 
Austrahan Research Council (Austraha), CNPq (Brazil) and CONICET (Argentina). 
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3.2. Deconvolution 

In order to enhance our sensitivity to the detailed structure of any extended emission, 
we removed the effects of the residual seeing halo from both the speckle SAA and the 
adaptive optics images by applying a deconvolution algorithm to selected regions containing 
each of the Northern Arm sources. After investigating a variety of deconvolution routines, 
including bhnd deconvolution (IDAC, Jefferies & Christou 1992; Christou et al. 1999), 
Lucy- Richardson (Lucy 1974; Richardson 1972), maximum-entropy and maximum-likelihood 
algorithms (Press et al. 1983), we concluded that the IDAC deconvolution routine provides 
the most robust results, based on the consistency of the observed extended morphology 
between maps taken in different years. This physically constrained algorithm^ was initially 
described by Jefferies & Christou (1992) and its application to AO imaging by Christou et al. 
(1999). In summary, IDAC is a multi-frame blind deconvolution package which makes use of 
multiple observations of the same source, where the differences between the observations are 
assumed to be due solely to differences in the point spread function (PSF), in order to solve 
for the underlying object distribution function as well as the PSF. The solution is driven 
by minimizing, using conjugate gradient techniques, an error metric that characterizes the 
difference between the observations and the convolution of the target and PSF estimates 
(See Appendix A for more details.). 



4. Results 
4.1. Morphology 

Figure 1 shows contour plots of the 0'.'5x0'.'5 region for the sources included in this 
study, all of which are clearly extended when compared to a neighboring PSF. IRS IW, 5, 
8, and lOW show signs of the horseshoe morphology expected from a bowshock. IRS 21 
shows an azimuthally symmetric extension, which is consistent with both the morphology 
and interpretation of a face-on bowshock as was deduced from the non-deconvolved images 
in Tanner et al. (2002). The bowshock hypothesis is, therefore, well supported based on a 
simple qualitative analysis of the images. 

To estimate the extent of the bowshock structures, the central compact component 
is subtracted from each deconvolved image (see Appendix and Figures 2 and 4). An initial 
estimate of the angular extent of the diffuse emission, rj, is then determined from the angular 



The current version of the code, implemented by E.K. Hege and M. Chesalka, is available from the 
Center for Adaptive Optics webpages at http://cfao.ucolick.org/software/idac/ 
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distance between the centroid position of the central source (prior to subtraction) and the 
position of the peak in the diffuse emission (after subtraction). As illustrated in Figure 2, 
this distance is used to quantify the two-dimensional structure by extracting radial profiles 
averaged over an azimuthal angle, ^a^,g(rad)=0.05/rj(arcsec). This 9avg corresponds to 25, 
25, 10 and 25 degrees for IRS IW, 5, 8 and lOW, respectively. For those radial profiles 
along which the peak is >30% of the overall maximum of the diffuse emission, the radius at 
which each profile reaches its peak intensity is recorded and its corresponding uncertainty is 
estimated by fitting a gaussian to the radial profile (see Table 2). 

The extent of these bowshocks allows us to apply models of the shape of the shock 
front to the observed morphologies. If the bowshocks are being produced from the flow of 
dust and gas along the Northern Arm encountering the wind from the central star, then a 
plane-parallel flow incident on a spherical wind would be the appropriate model. While the 
dust density of the Northern Arm is quite complex over large scales, we assume that over 
small scales it can be modeled as a constant density medium. In this case, the shape of 
a bowshock can be described analytically as R^-d^O) = RoCSc9-\/3{l — 9cot9), where i?o is 
the stand-off distance and 9 = cos~^9 ■ 9o is the angle between any direction 9 (unit vector) 
and the unit vector 9o describing the direction of the relative velocity between the central 
star and the incident medium, Vrei (Wilkin 1996; Canto et al. 1996). In order to fit our 
two-dimensional data set to this three-dimensional model, we consider only the points that 
lie in the plane defined by Vrei and Vrei x s, where s is the line of sight. In this plane, 
R2-D = {Rs-Dvrei^o, Ra-Dv^^ixdvreixs), whcrc thc i is a unit vector in the direction given 
by the subscript. Next, we project this two-dimensional model onto the plane of the sky, 
according to the inclination angle, i, the angle between Vrei and the line-of-sight. This gives 

R2 —D,sky 

{R^^Dv,„i9o,skySini, Rs ), where 9o,sky is 9o projected onto the plane 

of the sky ^. We use minimization to find the best fitting bowshock model to the data 
presented in Table 2. This provides estimates of the standoff distances, R^, and the projected 
orientation of 9o, the latter of which is described by the position angle on the plane of the 
sky, PA=tan"^(RQ,(^ = 9o)/Rs{9 = §0)) and the inclination. 

Figure 4 shows the best fitting bowshock models for IRS IW, 5, 8, and lOW while 
Table 3 gives the best-fitting parameters and their uncertainties, which are estimated from 
the change in the value of the parameter which results in a variation of the value by 1. 
When we try to fit this simple bowshock model to IRS IW, it was discovered that a better 



^This procedure projects onto the line of sight only the intersection of the bowshock with a plane passing 

through the central star, the apex of the bowshock and the vector s. It is thus only an approximation to 
the projected brightness distribution, a better approximation for which would be the curve represented by 
those lines of sight which are tangent to the 3-dimensional bow shock surface. 
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fit to the data could be achieved if the position of the apex of the model is displaced 200 
AU to the North. While this would reduce the value of the standoff distance by the same 
amount, the value of the inclination of the best fitting model remains 90 degrees. We plot 
both model fits in Figure 3 for comparison and use a value of 90 degrees for the inclination 
of the bowshock in subsequent calculations. The shift in the apex is of no concern as it could 
be explained as an offset in the position of the photocenter due to winds from neighboring 
sources while the poor fit to the data in the wings of the bowshock for the unshifted apex is 
due to the simplicity of the bowshock model which does not account for a density gradient 
or the magnetic field within the Northern Arm. 



4.2. Proper Motion of the Northern Arm Sources 

In order to investigate the relative motion occurring between the embedded sources and 
the Northern Arm, the proper motions of IRS 21, lOW, and IW are measured from the 
multi-epoch 2.2 /im speckle data set. The proper motions are estimated from a linear least- 
squares fit through their positions as a function of time. The uncertainties on the positions 
are estimated by assuming that the errors are equal for each source and in each direction and 
by requiring that their velocity fits have a equal to one; this method produces positional 
uncertainties for IRS 21 that are consistent with those found in Ghcz ct al. (1998), who 
explicitly measured the positional uncertainties. Table 4 gives the proper motions and their 
uncertainties as well as the positions of the sources at a specified epoch, while Figure 5 shows 
the positions of the sources with respect to time. Compared to our initial report of IRS 21's 
proper motion (Ghez et al. 1998), the value given here is consistent, but has a factor of 
five smaller uncertainty due to the additional epochs of data. For the other sources, these 
are the first reported estimates of their proper motions. Table 4 also provides the positions 
and proper motions of the Northern Arm sources estimated by Ott (2004). We will use the 
uncertainty- weighted average of these two sets of observations in the subsequent calculations 
in this paper. 



4.3. Model of the Northern Arm 

In order to explore the origin of the bowshocks as is done in §5, we need a kinematic 
model of the flow of the Northern Arm. Table 5 presents the values used, which come 
from the Paumard et al. (2004) three-dimensional flow model which incorporates a series of 
non-planar Keplcrian orbits. The velocities given in Table 5 represent the proper motion of 
the material at the position of each Northern Arm source as predicted by this model. The 
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uncertainties provided are the half range in each parameter which results in a 2(7 variation 
in the value. 

5. Discussion 

The observed extended emission demonstrates that the extreme redness of these sources 
is, in part, caused by dust that is not intrinsic to the source, but is rather from the postshock 
region of a bowshock, as suggested by Tanner et al. (2002). With the measured bowshock 
geometry and proper motions for some of the Northern Arm sources, as well as the model- 
derived velocity vectors of the Northern Arm material, it is now possible to investigate 
the origin of the bowshocks. A comparison of the measured proper motions of the sources 
(Table 4) with the flow of the Northern Arm (Table 5), shows that they do not agree. These 
observations, therefore, further rule out the hypothesis that the Northern Arm sources are 
recently formed stars that are being transported inward in the Northern Arm flow, thus 
relieving the challenges presented by the short dynamical timescales associated with the 
Northern Arm. The lack of any bowshock alignment toward the direction of Sgr A* allows 
us to rule out the formation of the bowshocks as being due to the wind produced by the black 
hole. We, therefore, conclude that the observed bowshocks are generated by an interaction 
with the Northern Arm flow, as opposed to any other medium. We now consider a scenario 
in which the embedded sources are interacting with this material and can calculate the 
relative velocity vectors between the sources and the flow. In this case, the PAs of the 
resultant velocity vectors should correspond to the predicted PAs for bowshocks caused by 
an interaction with the Northern Arm. These values are given for IRS IW and lOW in 
Table 6 and are consistent with the values based on the morphology of the bowshocks given 
in Table 3. While the observed PA of IRS lOW is marginally consistent with what is expected 
due to the interaction with the Northern Arm, it is extended in the direction of the known 
red giant, lOEE, suggesting a possible interaction with the winds from that star. 

Standoff distances measured for the Northern Arm sources provide some insight into 
the spectral type of the underlying star. The luminosity of these sources (10^~^ L©, Tanner 
et al. 2001; Becklin et al. 1978) suggests that the most likely source types for the central 
objects are AGB stars (v„,,=10 km s~^ and m=10~^ Mo/yr, Lamers & Cassinelli 1999), main 
sequence O stars (vtu=2000 km s^^ and m=10^^ MQ/yr, Prinja, Barlow and Howarth 1990; 
Howarth and Prinja 1989) or WR stars (v,„=1000 km s"^ and m=2xl0^^ Mq yr~^, Abbott 
& Conti 1987). Table 6 lists the standoff distances [Ragb^ Rq) and Rwr respectively) for all 
these stars based on the standoff distance relationship Ro — IHAxlQi^^m^^/'^Yw'^Y'^i^ix^Jj^n^^'^ 
cm (Weaver et al. 1977), where Ro is the standoff distance, is the stellar mass loss rate 
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in units of 10~^ M© yr""-*^, is the stellar wind velocity in units of 10^ cm s~^, Yrei* is the 
motion of the star relative to the Northern Arm in units of 10^ cm s~^, iih is the mean 
molecular weight, and is the gas density of the Northern Arm in units of cm^"^, where we 
use a value of 3x10^ cm~^ as derived by Tanner ct al. (2002). Only those standoff distances 
expected for the WR star scenario are consistent with the observed standoff distances. These 
large stand-off distances give the Northern Arm sources wind properties that are similar to 
the He I emission-hne stars observed in and around the central parsec (Krabbe et al. 1995; 
Genzel et al. 2000; Paumard et al. 2001, 2003). 

The distinct bowshock structure seen around IRS 8, combined with its projected distance 
from Sgr A*, has led to the hypothesis that it has been flung out of the central parsec by a 
gravitational interaction with the supermassive black hole (Geballe et al. 2004). Since the 
diffuse emission associated with IRS 8, which has been associated with the ionized gas at the 
northern tip of the Northern Arm (Lacy et al. 1979), is infalling toward Sgr A*, the large 
observed standoff distance between the star and the gas and dust can be explained by the 
star's large projected distance from Sgr A*. (40"). The large projected distance of IRS 8 
from the supermassive black hole means a smaller infall velocity for the surrounding diffuse 
material. As a result there is a smaller relative velocity between the star and diffuse material 
producing a larger standoff distance around IRS 8 than that seen around the embedded 
sources nearer to Sgr A*. This ehminates the requirement for the star itself to have a large 
proper motion, as is implied by the hypothesis that is has been flung out of the central parsec 
(Geballe et al. 2004). It has also been suggested that the near-infrared emission from the 
bowshock around IRS 8 is due to thermal shock emission from small grains (Geballe et al. 
2004). However, we have attributed the near-infrared emission around IRS 21 to scattered 
light, using a model of its SED (Tanner ct al. 2002). Both IRS 21 and IRS IW have large 2.2 
/im polarizations which are larger than that expected from the intervening ISM. They also 
have polarization position angles at 2.2 jim. which arc deviant from the position angles of 
surrounding stars (Ott et al. 1999). Both scattered and thermally emitted light from grains 
aligned in the asymmetric outflows of these stars provide an explanation for the observed 
near-infrared polarization properties of these stars. 

By modeling the asymmetric dust morphology of the Northern Arm sources as bow- 
shocks from stars with wind velocities on the order of 1000 km s~^, we can add IRS IW, 
low, and 21, and most probably IRS 5 and 8, to the population of hot, windy stars within 
the central parsec having winds and mass-loss rates resembhng those of WR stars. Within a 
24x24 arcsecond region around the central parsec, a total of 15 stars have been conflrmed as 
being candidate WR stars with strong winds (v^j„(i=300-1000 km s~^) through their broad 
He I emission-line profiles (Paumard et al. 2003, 2001; Figer, private comm.). With the 
width of the Northern Arm estimated from the FWHM of its emission at 12.5 /im (Tanner 
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et al. 2003) and the depth of the Northern Arm estimated from the half range of the hne-of- 
sight distances of IRS IW, 5, lOW and 21 taken from the model of Paumard et al. (2004), 
we estimate a volume filling factor of the Northern Arm of 10% within the 40x24 arcsecond 
field of view of the Paumard et al. (2001, 2004) He emission-line survey. Hence, given the 
known population of He I emission-line stars, we expect just a couple luminous stars to be 
coincident with the Northern Arm. Thus, the presence of three such stars (IRS 21, IW, 
and low) and probably more (IRS 2 and 5) within this same region is slightly larger than 
expected given the observed population of He I emission-line stars suggesting the population 
of windy stars has been underestimated to date. 

The presence of OB and WR stars, which are believed to be massive and quite young, 
within a parsec of the central supermassive black hole is difficult to explain. Strong tidal 
forces and low present gas densities make in situ star formation unlikely, while formation 
at larger radii and subsequent inward migration is challenging given the short dynamical 
timescales. Nonetheless, clear signatures of tangential anisotropy in the He I emission-line 
stars suggest that these are indeed a very young, unrelaxed population. Furthermore, for the 
subset of those stars with three-dimensional velocity vectors and two-dimensional position 
vectors, it has been shown that the orbits of the He I emission-line stars are consistent 
with at least one and possibly two rotating planes (Levin & Beloborodov 2003; Genzel et 
al. 2003). The identification of these planes, however, is model dependent since the He I 
emission-line stars have only five of the six position- velocity coordinates necessary to define 
the star's orbital planes. 

The Northern Arm sources offer the first opportunity to test the validity of these planes 
since all six phase-space values are known. Based on the proper motion data alone, IRS IW 
and low are found to be moving in a clockwise direction around Sgr A*. We estimate the 
radial velocities of IRS IW and lOW (see Table 4) from the derived inchnation angle of the 
bowshock, the velocity of the gas, and the relative proper motion of the gas and the star. The 
line-of-sight distance from Sgr A* is taken to be the same as the distance of the Northern Arm 
as given by the model from Paumard et al. (2003). Adding the measured proper motions 
and on-sky offsets from Sgr A*, we implement all six phase-space coordinates to estimate 
the angles between the planes described by orbits of IRS IW and lOW (see Table 7) and the 
planes of He I emission-fine stars given in Genzel et al. (2003). Using the same line-of-sight 
normalized angular momentum as Genzel et al. (2003), J^/J2(max), we determine that both 
IRS IW and lOW are on clockwise orbits. IRS IW or lOW both fie in the same plane as 
the clockwise plane of He I emission-fine stars (see C angle in Table 7) , suggesting that the 
Northern Arm sources, also being He I emission-line stars, share the same orbital plane as 
the known He I emission-line stars indicating that these disk hypotheses are valid for these 
sources. 
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6. Summciry and Conclusions 

With high-resolution images reveahng bowshocks around a population of luminous, em- 
bedded Northern Arm sources, we are able to conclude that they are all windy Wolf-Rayet 
stars and are thus likely to be members of the population of broad-line He I emission-line 
stars observed in the nearby IRS 16 cluster. They are not YSO's as previously speculated 
(Krabbe et al. 1995). The resolved bowshock structures are created from the interaction 
of the intrinsic stellar winds and the flow of material along the Northern Arm. While IRS 
IW and low show stellar wind properties and kinematics similar to the known population 
of clockwise orbiting He emission-line stars, their fully determined orbital parameters are 
coplanar with these sources. This supports the hypothesis that the He emission-line stars 
occupy at least one coherent orbital plane. The Northern Arm sources offer a unique op- 
portunity to test these plane models, so the addition of the known orbits of the remaining 
embedded sources (IRS 5, 8 and 2) would be beneficial. 

To further classify these Wolf-Rayet stars and attempt to directly measure their radial 
velocities, deeper spectroscopic observations with high spatial resolution in the near-infrared 
would be worthwhile, with the goal of finding the obscured spectroscopic signature of the 
stellar photosphere or wind. Chiar et al. (2003) observe the 6.2 /xm PAH (circumstellar 
amorphous carbon) feature toward the Quintuplet sources, which is known to be observed 
in intrinsic circumstellar dust shells around late type WR stars (DWCLs). If any absorp- 
tion lines intrinsic to the star are detected, the resulting radial velocity measurements of 
the Northern Arm sources would also be a powerful tool in understanding the dynamics of 
the He emission-line star population since their location within the Northern Arm allows us 
to estimate their full phase space coordinates with no assumptions about their line-of-sight 
distance. There have been recent 3.8 /xm observations of additional sources with bowshock 
morphologies within the central parsec, though these sources lack a 2.2 //m extended coun- 
terpart (Clenet et al. 2003). We are just beginning to understand the significance of these 
sources and with a more complete census of the region will be able to further investigate the 
population of windy stars within the central parsec through their interaction with the dust 
and gas within the mini-spiral. 

The work was supported through NSF grants No. AST-9988397 by the National Science 
Foundation Science and Technology Center for Adaptive Optics, managed by the University 
of California at Santa Cruz under cooperative agreement No. AST-9876783 and the Packard 
Foundation. We thank Thibaut Paumard for the use of his model of the fiow of the gas 
along the Northern Arm and also Eric Becklin for his many useful comments and ideas. 
Data presented herein were obtained at the W.M. Keck Observatory, which is operated 
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Cahfornia and the National Aeronautics and Space Administration. The Observatory was 
made possible by the generous financial support of the W.M. Keck Foundation. The authors 
wish to recognize and acknowledge the very significant cultural role and reverence that the 
summit of Mauna Kea has always had within the indigenous Hawaiian community. We are 
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A. IDAC and PSF Subtraction 

Since IDAC solves for both the object's spatial intensity distribution function as well as 
the PSF, it has more flexibility in the way it can be run than most deconvolution algorithms. 
Prior to deconvolution, an initial estimate of the target image is obtained from the mean 
of the input images, and an initial estimate of the PSF is created from a two-component 
Gaussian fit to a reference star in the field. The latter is a reasonable model of the PSF, 
modeling both the central diffraction-limited core and the broader, lower-intensity residual 
seeing halo. The IDAC algorithm is initially run on the data cube of three to five images 
for ~200 iterations, allowing both the target and PSF images to vary. During this process, 
a point source in the field does not reduce to a delta function in the target reconstruction 
image, requiring that the reconstructed PSF be reconvolved with an isolated point source 
in the target image to produce a new estimate for the PSF. The algorithm is then restarted 
for a couple hundred more iterations with the new PSF and target image from the previous 
deconvolution. This procedure is repeated a couple of times until a good PSF is obtained for 
each frame within the data cube. The algorithm is then run for about 20 iterations with the 
PSF held fixed to remove any residual PSF structure from the target estimate. Finally, the 
algorithm is run for a few more iterations with both the target and PSF permitted to vary 
in order to relax the solution thus producing the final deconvolved image and PSF estimate. 

IDAC has a known tendency to produce point sources whose size depends very modestly 
on their brightness. Over five magnitudes, the core size is observed to increase by 35% 
(Christou et al. 2004 in prep). This effect is negligible compared to the stand-off distances 
observed in the Northern Arm sources, but needs to be accounted for in order to carry out 
PSF subtraction (§4.1). Since the Northern Arm sources are all brighter than their nearby 
point sources, the later are convolved with a gaussian to provide the best match for PSF 
subtraction. In order to avoid contamination from the bowshock emission the comparison 
of core sizes is carried only out to source radii that correspond to 70% of the peak intensity 
level. Figure 4 shows images of the Northern Arm sources deconvolved with IDAC and PSF 
subtracted using the Gaussian convolved point source. 
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Fig. 1. — Contour 2.2 /im images of the 0'/5x0'.'5 region around the Northern Arm sources. 
Column 1 shows the raw images, Column 2 is a corresponding PSF taken from the same 
image, Column 3 shows the images deconvolved using the IDAC deconvolution algorithm 
(Jcfferies & Christou 1992). The contours plotted represent 90-20% of the peak value in 
10% intervals. North is up and East is to the left. 

Fig. 2. — a) Contour image of IRS 5 along with lines depicting the arcs over which radial 
profiles were averaged. North is up and East is to the left, b) Plots of those radial profiles 
for which the peak is within 70% of the maximum value of all the profiles. 

Fig. 3. — Contour images of the two best fits of the bowshock model to the radii extracted 
from the image of IRS IW. Left: The best fit derived after shifting the apex of the bowshock 
200 AU to the North. Right: The best fit derived without shifting the apex. 

Fig. 4. — Center: Greyscale image of the central 20x20 arcseconds of the galactic center at 
12.5 ixni presented in Tanner et al. (2002) overplotted with the symbols representing the 
position of the Northern Arm stars (+), the broad-hne He emission-hne stars (diamonds) 
and narrow-hne He emission-hne stars (*) from Paumard et al. (2001, 2003). Sides: Contour 
images of IRS IW, 5, 8 and lOW from the deconvolved Keck/NIRC and Gemini/AO 2.2 pim 
data after PSF subtraction. Overplotted as crosses in the contour images are the positions of 
the peaks in the radial profile used to find the best fitting bowshock models (superimposed 
curve). North is up and East is to the left. 

Fig. 5. — Plots of the a and S position of IRS 21, IW, and lOW, as a function of time, along 
with the least squares fit to the data. 
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Table 1. Summary of Observations 



lelescope / Camera 


Date 


Sources 


Kork T/NTRC 




21" 


Keck I/NIRC 


1996 Time 


21" 


Keck I/NIRC 


1 997 Mav 


21" 


Keck I/NIRC 


1998 May 14-15 


1W,10W,21" 


Keck I/NIRC 


1998 August 4-6 


1W,10W,21" 


Keck I/NIRC 


1999 May 2-4 


1W,10W,21" 


Keck I/NIRC 


1999 July 24-25 


1W,10W,21" 


Keck I/NIRC 


2000 May 19-20 


1W,10W,21" 


Keck I/NIRC 


2000 July 19-20 


1W,10W,21" 


Keck I/NIRC 


2001 May 9 


1W,10W,21" 


Keck I/NIRC 


2001 July 28-29 


1W,10W,5,21" 


Keck I/NIRC 


2002 April 23-24 


21 


Keck I/NIRC 


2002 May 23-24 


21 


Keck I/NIRC 


2002 July 19-20 


21 


Keck I/NIRC 


2003 April 21-23 


21 


Keck I/NIRC 


2003 July 22 


21 


Keck I/NIRC 


2003 September 7-8 


21 


Gemini/Hokupa'a-QUIRC 


2000 July 6 


8 



^2.2 jjLTo. observation of IRS 21 were presented originally by Tanner et al. 2002. Deconvo- 
lutions of these observations are introduced here along with proper motion measurements. 
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Table 2. Azimuthal Dependence of Bowshock Morphologies Around the Northern Arm 

Sources 



Source 


PSF Source" 


PA 


Intensity 


Radius [AU] 


PA 


Intensity 


Radius [AU] 


IRS IW 


1NE(3) 


10 


0.530+0.005 


470+80 


150 


0.591+ 0.003 


710+140 






30 


0.556+0.003 


460+90 


250 


0.450+0.003 


720+90 






50 


0.520+0.003 


470+100 


270 


0.762+0.003 


620+80 






70 


0.521+0.003 


470+100 


290 


0.520+0.003 


580+80 






90 


0.749+0.003 


490+90 


310 


1.000+0.003 


560+80 






110 


0.768+0.003 


460+100 


330 


0.969+0.003 


500+80 






130 


0.610+0.003 


570+120 


350 


0.774+0.004 


480+80 


IRS 5 


-1.41+0.21 


-10 


0.856+ 0.004 


560+100 


90 


0.654+ 0.004 


510+90 






20 


0.637+ 0.004 


470+100 


115 


0.840+0.004 


520+90 






40 


0.351+ 0.004 


440+100 


140 


1.000+0.004 


610+90 






70 


0.404+ 0.004 


460+90 








IRS 8 


+1.52-5.0 


35 


0.437+ 0.004 


2030+200 


270 


0.408+ 0.002 


2070+140 






40 


0.629+ 0.003 


1980+180 


280 


0.492+0.003 


1940+160 






50 


0.889+ 0.003 


1940+160 


290 


0.435+0.003 


1820+180 






60 


1.000+ 0.003 


1920+160 


300 


0.426+0.003 


1700+160 






65 


0.946+ 0.003 


1910+160 


310 


0.441+0.003 


1600+160 






75 


0.818+ 0.003 


1910+160 


315 


0.428+0.003 


1520+160 






80 


0.656+ 0.003 


1960+160 


320 


0.412+0.003 


1500+160 


IRS low 


lOEE 


30 


0.526+ 0.004 


460+100 


120 


1.000+0.004 


410+100 






60 


0.675+ 0.004 


400+80 


150 


0.883+0.004 


520+120 






90 


0.965+ 0.003 


430+190 


180 


0.628+0.005 


760+50 



^PSF Sources for IRS 5 and 8 are labeled as the North-East offsets from the extended 
sources in arcseconds. 
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Table 3. Measured Bowshock Properties 



Source 


Standoff Distance [AU] 


Inclination [°] 


PA [°] 


IRS IW" 


210±140 


90±10 


30±20 


IRS IW 


460±160 


90±10 


30±30 


IRS 5 


450±60 


80/100±20 


60±20 


IRS 8 


1380±300 


90±10 


0±10 


IRS low 


420±60 


70/110±20 


85±20 


IRS 21" 


670±40 







^Values derived from the best fitting bowshoclc model 
after moving the apex 200 AU to the North 

^IRS 21 standoff distance as estimated by Tanner et al. 
(2002) 
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Table 4. Kinematics of the Northern Arm Sources 



Source 


to 


ato 




z 


a 

a 




z 




year 


arcseconds 


arcseconds 


arcseconds 


km s~^ 


km s~^ 


km s~^ 


IRS IW 


1999.164 


5.283±0.007 


0.571±0.008 


7.537±1.536 


-170±30 


320±60 


20±110 


IRS low 


1999.905 


6.492±0.009 


5.108±0.011 


6.287±1.447 


-160±110 


10±210 


175±130 


IRS 21 


2000.207 


2.383±0.006 


-2.713±0.007 


4.455±1.080 


-300±30 


100±30 




IRS IW 




5.251±0.007 


0.580±0.009 




-100±29 


216±37 


20±90 


IRS low' 




6.507±0.001 


5.117±0.001 




51±15 


176±15 


215±160 



'^Position and velocity on the plane of the sky measured from 2.2 /xm images. Additional 
systematic sources of error in the positions arise from uncertainties in the plate scale and 
image rotation. 

is estimated from the models of the Northern Arm with the assumption that these 
sources are embedded in the Northern Arm. In this convention, the line of sight axis increases 
away from the observer. 

is estimated from the bowshock modehng described in § 4.1. 

'^Position and velocity on the plane of the sky from Ott (2004). We use the weighted 
average of these two measurements in subsequent calculations. 



Table 5. Properties of the Flow of the Northern Arm 



Location 


Va 


V5 






km s~^ 


km s~^ 


km s"-*" 


IRS IW 


-70±50 


-210±70 


20±50 


IRS 5 


10±130 


-140±190 


110±60 


IRS low 


-20±70 


-190±100 


80±60 


IRS 21 


-130±60 


-320±50 


-90±20 
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Table 6. Predicted PAs and Standoff Distances for the 
Embedded Sources 



Source 




^AGB 








degrees 


AU 


AU 


AU 


IRS IW 


-10±10 


40±10 


60±10 


560±40 


IRS low 


10±10 


50±10 


70±10 


680±30 


IRS 21 


-20±10 


40±10 


55±10 


550±30 



'^This is the measured PA of the relative velocity between 
the sources and the Northern Arm flow on the plane of 
the sky and is expected to correspond to the PA of the 
bowshock if the bowshock is created by an interaction of 
the star's wind and the Northern Arm material. 

''The uncertainties only incorporate our measured un- 
certainties for the relative velocities and stellar proper mo- 
tions. 



Table 7. Estimated Orbital Planes 



Source 


Observed i 
degrees 


Orbit i 
degrees 


Orbit Q 

degrees 


Jz/ Jz{'mcix) 


C angle" 
degrees 


CC Angle" 
degrees 


IW/lOW an! 
degrees 


IRS IW 


90 


123±6 


149±14 


0.93±0.16 


27±35 


162±60 


47±53 


IRS low 


70 


129±27 


90±44 


0.61±0.10 


25±114 


131±39 




IRS low 


110 


121±15 


-155±38 




73±31 


139±43 


47±43 



C and CC stand for clockwise and counter clockwise. 
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